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Dynamical models of prototype gravastars were constructed in order to study 
their stability. The models are the Visser- Wiltshire three-layer gravastars, in which 
an infinitely thin spherical shell of stiff fluid divides the whole spacetime into two 
regions, where the internal region is de Sitter, and the external is Schwarzschild. It 
is found that in some cases the models represent the "bounded excursion" stable 
gravastars, where the thin shell is oscillating between two finite radii, while in other 
cases they collapse until the formation of black holes. In the phase space, the region 
for the "bounded excursion" gravastars is very small in comparison to that of black 
holes, but not empty. Therefore, although the existence of gravastars cannot be 
excluded from such dynamical models, our results do indicate that, even if gravastars 
indeed exist, they do not exclude the existence of black holes. 
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I. INTRODUCTION 

Black holes are well known and accepted objects not only on the scientific community 
but also on the general public. This is, for example, attested to partially by the phenomenal 
success of the well-known film, Black Holes (New River Media, 1998), directed by Pappi 
Corsicato, and by several best selling popular books, such as Kip S. Thorne, Black Holes and 
Time Warps: Einstein's Outrageous Legacy (W.W. Norton & Company, 1995), and Stephen 
Hawking, A Brief History of Time (Bantam Books publisher, 10th edition, 1998). However, 
the real detection of black holes highly demands their detailed explorations, including the 
forms of gravitational waves emitted by them. In principle, such a detection can never be 



conclusive and is fundamentally impossible [1], although there are many strong astronomical 
evidences for their existence. 

The simplest example of black holes is the spherically symmetric vacuum solution of the 
Einstein field equations, 

ds 2 = -f(r)dt 2 + -^--dr 2 + r 2 (d6 2 + sin 2 9d<p 2 ) , (1.1) 
f{r) 



where the function f(r) is given by 



/M = l-^, (1.2) 
r 



in the units where c = 1 = G. The metric is singular at both r = and r = r g = 2M. 
However, the nature of these singularities is different. In particular, the one at r = is 
generic and the spacetime curvature diverges there, while the one at r = r g is a coordinate 
one and can be made disappear after proper coordinate transformations. A classical point 
test particle will freely fall through r = r g without experiencing anything special. Such a 
surface is usually called as an event horizon. The studies of the properties of such horizons 
are fundamentally important, and it is found that there may exist a deep connection between 
gravity and thermodynamics The discovery of the quantum Hawking radiation [3] and 
ihe black hole entropy which is proportional to the area of the event horizon of the black hole 
4| further supports this idea. This interesting relation was first manifested when Jacobson 
derived the Einstein field equations from the first law of thermodynamics (FLT) by assuming 
the proportionality of the entropy and the horizon area for all local acceleration horizons jsj. 
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For static spherically symmetric and stationary axisymmetric space-times, Padmanabhan 



et al showed that Einstein's equations at the horizon give rise to the FLT |6|, |7|. Such 



BEL 

considerations were further generalized to the Lovelock gravity |8j]. In [9j], on the other 
hand, the gravitational field equations (GFE's) for the nonlinear f(R) theory were derived 
from the FLT by adding some non-equilibrium corrections. 

For the de Sitter space there also exist Hawking temperature and entropy associated with 
the cosmological event horizons and its thermodynamic laws [loj]. In this space, the event 
horizons coincides with the apparent horizon. For more general cosmological models, event 
horizons may not exist, but apparent horizons always do, so it is possible to have Hawking 
temperature and entropy associated with apparent horizons. Along this line, the connection 
between the FLT of apparent horizons and the Friedmann equations in Einstein's th eory 



with/without the Gauss-Bonnet term, as well as in the Lovelock theory, were found 11]. 
More recently, it was found that this is also true for the braneworld cosmology 12]. In 13), 
with the help of a new mass-like function, it was shown that the FLT of apparent horizons 
in equilibrium can be derived from the Friedmann equations in various theories of gravity, 
including the Einstein, scalar-tensor, nonlinear f(R), and Lovelock. 

Despite of all these theoretical and observational successes, a number of paradoxical 
problems reparging to black holes also exist [jjj , which frequently motivate authors to look 
for other alternatives, in which the endpoints of gravitational collapse are massive stars 



without horizons. Example of such models include gravastars [l5|], Bose superfluid 16], 
and black stars [l^j], to name only few of them. Among these models, gravastars have 
received particular attention recently [181 ]. partially due to the tight connection between 
the cosmological constant and a currently accelerating universe [19j, although very strict 
observational constraints on the existence of such stars may exist 20j. In the original 
model of Mazur and Mottola (MM) 15], gravastars consist of five layers: an internal core 
< r < r%, described by the de Sitter universe, an intermediate thin layer of stiff fluid 
ri < r < V2, an external region r > 7*2, described by the Schwarzschild solution (II. ip - 
(II. 2j) . and two infinitely thin shells, appearing, respectively, on the hypersurfaces r = n and 
r = 7"2- By properly choosing the free parameters involved, one can show that the two shells 
can have only tensions but with opposite signs 15]. Visser and Wiltshire (VW) argued that 



such five-layer models can be simplified to three-layer ones 



21(, in which the two infinitely 



thin shells and the intermediate region are replaced by one infinitely thin shell, so that the 



r = 



(a) 



r = a(x) 



(b) 



FIG. 1: (a) The Penrose diagram for the Schwarzschild vacuum solution. Region II is the region 
where r > a(r). (b) The Penrose diagram for the de Sitter solution. Region / is the region 
where r < a(r). An infinitely thin shell located at r = a(r) connects these two regions to form a 
dynamical spacetime of a prototype gravastar, as shown in Fig. [2j 

function f(r) in the metric (11.11) is given by 



where r = a(r) is a timelike hypersurface, at which the infinitely thin shell is located, and 



where an overdot denotes the derivative with respect to the proper time r of the thin shell. 
Therefore, in the region r > a(r) the spacetime is locally Schwarzschild, while in the region 
r < a(r) it is locally de Sitter. These two different regions are connected through a dynamical 
infinitely thin shell located at r = a(r) to form a new spacetime of gravastar [cf. Figs. Q] 
and % 

To study the dynamics of Eq. fll.4p . one can follow two different approaches: one is to 
prescribe a potential V(a) and leave the equation of state of the shell as derived, and the 
other is to prescribe an equation of state of the shell and leave the potential V(a) as derived. 
VW followed the first approach, and studied in details the case where 




(1.3) 



r denotes the proper time of the thin shell. The constant I = a/3/A denotes the de Sitter 
radius. On the hypersurface r = o(r) Israel junction conditions yield 




V(a ) = 0, V'(a ) = 0, V" (a Q ) > 0, 



(1.5) 
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FIG. 2: The Penrose diagram for the solution given by Eq, (|1.3p . The infinitely thin shell at r = a(r) 
connects the two regions I and // shown in Fig. Q]to form a spacetime of a prototype gravastar. 
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where a prime denotes the ordinary differentiation with respect to the indicated argument. 
If and only if there exists such an a for which the above conditions are satisfied, the model 
is said to be stable. Among other things, VW found that there are many equations of state 
for which the gravastar configurations are stable, while others are not 2l|. Carter studied 
the same problem and found new equations of state for which the gravastar is stable 22], 



while De Benedictis et al 



23| and Chirenti and Rezzolla 24| investigated the stability of the 



original model of Mazur and Mottola agaist axial-perturbations, and found that gravastars 
are stable to these perturbations. Chirenti and Rezzolla also showed that their quasi-normal 
modes differ from those of a black hole of the same mass, and thus can be used to discern a 
gravastar from a black hole. 

As VW noticed, there is a less stringent notion of stability, the so-called "bounded ex- 
cursion" models, in which there exist two radii a± and a 2 such that 



V( ai ) = 0, V'( ai )<0, V(a 2 )=0, V'(a 2 )>0, 



;i.6) 



with V(a) < for a E a 2 ), where a 2 > a±. 

In this paper, our purpose is twofold: (a) First, we construct three-layer VW dynamical 
models, some of which represent "bounded excursion" stable gravastars, and some represent 
the collapse of a prototype gravastar, where the final fate of the collapse is the formation of 
a black hole. (2) Second, in the phase space we compare the region of such stable gravastars 
with the one of black holes, and show explicitly that both of them are non-zero, although 
the former is much smaller than the latter. The rest of the paper is organized as follows: 
In Sec. II we shall study various cases, in which all the possibilities of forming black holes, 
gravastars, de Sitter, and Minkowski spacetime exist. In Sec. Ill we present our main 
conclusions. 



Before turn to the next section, we note some relevant work. In particular, t 
tional collapse of dark energy in the background of dark matter was studied in 



re gravita- 



251 ] . and it 



was found that when only dark energy is present, black holes are never formed. When both 



of them are present, black holes can 



De formed, due to the condensation of the dark matter. 



Similar results were obtained in 



26 



281 ] . Recently, such studies were further generalized to 



a homogeneous and isotropic expanding Friedmann-Robertson- Walker universe dominated 



by dark energy 



29]. 
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II. FORMATION OF GRAVASTARS FROM GRAVITATIONAL COLLAPSE 

To keep the ideas of MM as much as possible, we consider the thin shell as consisting 
of a stiff fluid, a = where a and i? denote, respectively, the surface energy density and 



tension of the shell. Then, we find that 21] 



" = ' °(") 4 ' (") 

where gq and are integration constants, and have dimensions of surface energy density 
and length, respectively. It can be shown that the potential appearing in Eq.f ll.4p now can 
be cast in the form, 

- yw If m 1 9 „ 4 R 2 mR 7 R 10 \ , . 

V(R) = - 2 (-l + - + w + rf# + — -— + -), (2.2) 

where 

m = 77173' R =kUz> ^ = 77773 ' ( 2 - 3 ) 
with k = 4vTaQO"o. Therefore, for any given constants m and L, Eq. fll.4p uniquely determines 
the collapse of the prototype gravastar. Depending on the initial value Rq, the collapse can 
form either a black hole, or gravastar, or a Minkowski, or a de Sitter space. In the last case, 
the thin shell first collapses to a finite non-zero minimal radius and then expends to infinity. 
To guarantee that initially the spacetime does not have any kind of horizons, cosmological 
or event, we must restrict Ro to the range, 

2m<R < L, (2.4) 

correspondingly r G (r g ,l). When m = = A, the thin shell disappears, and the whole 
spacetime is Minkowski. So, in the following we shall not consider this case, and begin with 
the one m = and A / 0. 

A. m = and A / 

In this case, the spacetime outside the thin shell is flat, and the mass of the shell com- 
pletely screens the mass of the internal de Sitter spacetime. From Eq. (12.21) we find that 

i / -1 p2 R 10 \ 
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FIG. 3: The potential V(R) for m = 0. The top line is for L = 2.0 > L c ~ 0.9588, the middle line 



is for L = L c , and the bottom line is for L = 0.5 < L c . 



Then, it can be shown that the equations V(R) = and V'(R) = have the explicit solution, 



5 /^ 8/3 ' 



3 V5 



1/2 



0.9588, R = R min = ( - ) . (2.6) 



1/3 



For L < L c the potential V(R) is strictly negative as shown in Fig. [3j As a result, if the star 
starts to collapse at R = Rq, it will collapse continuously until R = 0, whereby a Minkowski 
spacetime is formed, as shown by the bottom line in Fig. HI When L = L c ~ 0.9588, since 
Rq < L c , we can see that, similar to the last case, the star will collapse until the center 
R = and turns the whole spacetime into Minkowski, as shown by the middle line in Fig. HI 
For L > L c , the potential V(R) is positive between Ri and R2, where R\^ are the two real 
roots of the equation V(R,L > L c ) = with R 2 > R\ > 0. In this case, if the star starts 
to collapse with Rq < Ri, as can be seen from Fig. [3j it will collapse to R = 0, whereby 
a Minkowski spacetime is finally formed. If it starts to collapse with Rq > R2, it will first 
collapse to R = R2 and then starts to expand until R = 00, and the whole spacetime is 
finally de Sitter, as shown by the top line in Fig. HI 



B. A = and m / 

In this case, Eq. fj2.2H yields, 

. . 1 / 771 1 r, A . 

v w = 2{ 1 -g-m-" f *)' (2 ' 7) 



9 




FIG. 4: The motion of the shell R(t) vs the proper time, r, of the shell for m = 0. The top line is 
for L = 2.0 > L c ~ 0.9588, the middle line is for L = L c , and the bottom line is for L = 0.5 < L c . 

from which we find that the equations V(R) = and V'(R) = have the explicit solution, 

m = ™ c =\ Gf - ">-*-■ - (if - (£f- («) 

For m > m c the potential V(R) is strictly negative as shown in Fig. El Then, the collapse 
always forms black holes, as shown clearly by the bottom line in Fig. |6j For m = m c , there 
are two different possibilities, depending on the choice of the initial radius Rq. In particular, 
if the star begins to collapse with Rq > R c , the collapse will asymptotically approach the 
minimal radius R c . Once it collapses to this point, the shell will stop collapsing and remains 
there for ever, as can be seen from the middle line in Fig. [6j However, in this case this point 
is unstable, and any small perturbations will lead the star either to expand for ever and leave 
behand a flat spacetime, or to collapse until R — 0, whereby a Schwarzschild black hole is 
finally formed. On the other hand, if the star begins to collapse with 2m c < Rq < R c , as 
shown by Fig. El the star will collapse until a black hole is formed. For m < m c , the potential 
V(R) has a positive maximal, and the equation V(R, m < m c ) = has two positive roots 
Ri t 2 with R 2 > R± > 0. As in the last case, now there are also two possibilities, depending 
on the choice of the initial radius Rq. If Rq > R 2 , the star will first collapse to its minimal 
radius R = R 2 and then expand to infinity, whereby a Minkowski spacetime is finally formed, 
as shown by the top line in Fig. EJ If 2m < Rq < Ri, the star will collapse continuously 
until R = 0, and a black hole will be finally formed. 
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FIG. 5: The potential V(R) for A = (or L = oo). The top line is for m < m c ~ 10~ a286 , the 
middle line is for m = m c , and the bottom line is for m > m c . 
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FIG. 6: The motion of the shell R(t) vs its proper time r for A = 0. The top line is for m < m c ~ 
1CP ' 286 , the middle line is for m = m c , and the bottom line is for m > m c . 

C. rn^O and A / 



In this case, from Eq.f l2.2p we find that the equations V(R) = and V'(R) = have the 
solution of the form, m = m c (L) for any given L. The exact dependence of m c on L cannot 
be given explicitly. Instead, in the following we consider some representative cases. 
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FIG. 7: The potential V(R) for m = 10 a286 and L = 2.5 x 10 5 , which is always negative. 

1. m = lO" - 286 , L = 2.5 x 10 5 

If we set m = m c ~ io -0 - 286 and graduately turn on the cosmological constant, we find 
that the potential V(R) becomes completely negative, as shown by Fig. [71 Then, for any 
given i?o with 2m < Rq < L, the star will always collapse to form a black hole. Comparing 
it with the case m = m c and A = 0, we find that the presence of the cosmological constant 
makes the collapse more like to form black holes than gravastars, or any of the others. 



2. m = 1(T 4 , L = 2.5 x 10 5 

If we keep L fixed, i.e., L = 2.5 x 10 5 , and tun m downward, we find that, for m = 10 -4 , 
the potential takes the shape given by Fig. [HI from which we can see that V(R) = now 
has four real roots, say, Ri, where Ri+i > Ri- If we choose Rq > R4, then the star will 
first collapse to R = R4, and then expand to infinity, whereby a de Sitter space is finally 
formed. However, if we choose R2 < Ro < R3, the collapse will bounce back and forth 
between R = R2 and R = R3. Such a possibility is shown in Fig. [91 This is exactly 



the so-called "bounded excursion" model mentioned in [21], but was not studied there or 
somewhere else. Of course, in a realistic situation, the star will emit both gravitational waves 
and particles, and the potential shall be self- adjusted to produce a minimum at R — Rstatic 
where V (R = Rstatic) = = V (R = Rstatic) [cf. Fig. [TO] , whereby a gravastar is finally 
formed 21]. 
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FIG. 9: The development of R(t) vs the proper time r for m = 10 4 and L = 2.5 x 10 5 . 

Although such a shape of potential is very difficult to find in the phase space of m and L, 
its measurement is not zero. So, one cannot completely exclude the existence of gravastars. 
However, we do find that it is easy to find potentials that lead to the formation of black 
holes. Figs. fTTI and fT2l are for m = 1.0 and L = 3.0, and m = 10 and L = 3.0, respectively. 
From these figures we can see that, by properly choosing the initial radius Rq, the collapse 
always forms black holes. In contrary, in these cases gravastars cannot be formed no matter 
how to tun Rq. On the other hand, we are not able to find values of m and L for which only 
gravastars can be formed. 



A V(R) 
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FIG. 10: The potential V(R) for the formation of a gravastar in a realistic collapse, after the star 
settles down to the minimus point R = Rstatic, where V (R = Rstatic) = = V (R = Rstatic)- 




FIG. 11: The potential V(R) for m = 1.0 and L = 3.0. 
III. CONCLUSIONS 



In this paper, we have studied the problem of the stability of gravastars by constructing 
dynamical three-layer models of VW [21|, which consists of an internal de Sitter space, a 
dynamical infinitely thin shell of stiff fluid, and an external Schwarzschild spacetime. We 
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FIG. 12: The potential V(R) for m = 10" 4 and L = 3.0. 

have shown explicitly that the final output can be a black hole, a "bounded excursion" 
stable gravastar, a Minkowski, or a de Sitter spacetime, depending on the total mass m 
of the system, the cosmological constant A, and the initial position R of the dynamical 
shell. All these possibilities have non-zero measurements in the phase space of m, A and 
i?o, although the region of gravastars is very small im comparing with that of black holes. 
Therefore, even though the existence of gravastars cannot be completely excluded in these 
dynamical models, our results do indicate that, even if gravastars exist, they do not exclude 
the existence of black holes. 
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